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The normalized multiplicity moments and their relation 
with soft and hard processes in relativistic heavy ion collisions 
are analyzed in a general two-component model. It is found 
that, the strong fluctnations in binary collision number Nc 
in minimum-bias events can enhance the hard component, 
especially for the higher order moments. This enhancement 
can not be effectively described by modifying the participant 
number in the one-component model. 


I. INTRODUCTION 

The relativistic heavy ion collisions at SPS and RHIC 
may be the only way to create the extreme conditions 
necessary to produce a new state of matter — Quark- 
Gluon Plasma (QGP) in the laboratory ||,^. One can 
attempt to understand the energy density achieved in 
the collisions by studying the multiplicity and transverse 
energy distributions through hydrodynamic models [^. 
At SPS energies, the global quantities like average multi¬ 
plicity, multiplicity distribution and rapidity distribution 
can be well described by soft processes only, namely by 
the number of participant nucleons only Q|. However, 
at RHIG energies, the measured pseudorapidity density 
normalized per participant pair for central Au-Au colli¬ 
sions shows that 70% more particles are produced than 
at SPS H,^. This indicates that the yield of particles 
created by hard scattering processes becomes important 
at RHIG [Q,Q. One can decompose the multiplicity at 
fixed impact parameter into a soft component and a hard 
component as 

n = aNp + bNc , (1) 

where Np and Nc are the participant number and binary 
collision number, respectively. 

However, the two-component expression (|^ can be ef¬ 
fectively described by a simple power-law form, 

n = , a>l , (2) 

which is then similar to that measured at SPS j|]. A nat¬ 
ural question then arises: Gan one find other global ob¬ 
servables which are more sensitive to the hard processes 
than the multiplicity itself, and which can not be effec¬ 
tively described in the models with only soft processes? 

As is well known, the multiplicity moments are im¬ 
portant characteristics in multiparticle production. The 
properties of the multiplicity distribution can be com¬ 
pletely described by the normalized moments 


z - 2 3 . 


( 3 ) 


In Ref. Ci were investigated at SPS energies with 
a general wounded nucleon model. It was found that 
the normalized multiplicity moments are independent of 
the concrete behavior of elementary nucleon-nucleon col¬ 
lisions, but dominated by the normalized participant mo¬ 
ments 


Ci — Cip — 
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provided that the colliding nuclei are not too light. 

In this paper we investigate how the hard processes 
change the normalized multiplicity moments. We extend 
the study in Ref. to including the hard component. 
We will focus on the sensitivity of Q to the colliding 
energy, nuclear geometry and especially to the geometry 
fluctuations. 


II. MULTIPLICITY MOMENTS 


At fixed impact parameter &, the nuclear geometry of 
soft and hard processes is expressed in terms of Np and 
Nc 1^, respectively. 


Np{b) = J TAis) (l - 

-krB(b-s) (^1 - 

Ncib) = J (fs aNTAis)TB{h - s) , 


( 5 ) 


where is the nucleon-nucleon inelastic cross section, 
and Ta{s) and TB(b — s) are the local participant densi¬ 
ties in the plane orthogonal to the collision axis defined 
as 


Ta{s) = J dzpA{s,z) , 


TB(b-s) 


J dzpsih - s,z). 


( 6 ) 


If the average multiplicity distribution of each soft 
source is gp{np), and the average multiplicity distribu¬ 
tion of each hard source is gduc), the multiplicity distri¬ 
bution of an AB collision at impact parameter b is the 
supposition of the contributions of Np soft sources and 
Nc hard sources: 
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It will be seen later that the results in this paper are not 
concerned with the concrete form of gp(np) and gdnc)- 
Taking into account all the processes with different b, the 
final-state multiplicity distribution of the nucleus-nucleus 
collisions is 


P{n)= Y. p{Np,N,)GN,,NSn) . ( 8 ) 

Np,N^ 

where p{Np, Nc) is the distribution function of Np and 
Nc- Since Nc calculated with is a monotonous func¬ 
tion of Np, p{Np, Nc{Np)) is just the distribution function 
p{Np) which can be obtained from (||) as 


we derive the multiplicity moments (n*) in terms of the 
elementary soft and hard moments (np and (np and the 
nuclear geometry moments (-/Vp, (iV*) and {N^N^), 

in) = {Np){np) + {Nc){nc) , 

(np = {{N^} - (Np)) {upf + {Np){nl) 

+ ((TVP - {Nc)) (nc)^ + {Nc){nl) 

+ 2{NpNc){np){nc) , 

(nP = ((7VpP-3(iVp+2(iVp) (n^p 

+ 3 ((iVp - {Np)) {np){nl) + {Np){nl) 

+ {{Nl)-^Nl)+2{Nc)) {nc)^ 

+ 3 ((iVp - (iVp) (nP(np + (7VP(nP 
+ 3 {{N^Nc) - {NpNc)) {np)^{nc) 

+ 3 {{NpN^) - {NpNc)) {np){nc)^ 

+ i{NpNc){nl){nc) + i{NpNc){np){nl) , 
. ( 12 ) 


with the definition of the moments, 


p{Np, Nc{Np)) = p{Np) ^ b{Np)^^ . (9) 

We introduce now generating functions [ p^Jl^ 
F{6), fp{0) and /c(^*) for the whole system and the el¬ 
ementary soft and hard sources, 

F{6) = ^0"P(n) 

n 

= ^p(ivp[/p0)]^M/c(0)]^'= , 

Afp 

fp{e) = Y<^'^^9p{np) , 

Up 

fc{9) = Y ^”‘=ffc(nP , -1 < 0 < 1 . (10) 

ric 

Differentiating (0) with respect to 9 and making use of 
the relations, 

P(0)P=l=/P0)P=l=/c(0)k=l = l, 

^Qn0)\e^i = {n) , 

^/p(^)|0=i = (''^p)) 

^/c(0)|e=i = (w) , 

^P(0)|e=i = (n(n- 1)) , 

^2 

g^/p(^)l®=i = {np{np — 1)) , 

o2 

j^fc{9)\0=i = {nc{nc-l)) ,■■■ ( 11 ) 
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«) = E ^cff(«c) , 
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{n;) = Y Kp(^p) - 

Np>Nmirt 

(iVp = ^ N),{Np)p{Np) , 

Np>Nmir, 

{n;ni)= Y N;N^ciNp)piNp) , (13) 

Np>Nmin 


where we have used the minimum participant number 
Nmin to select events. Nmin — 2 means minimum-bias 
events and very large Nmin corresponds to central events. 

With the known multiplicity moments, the normalized 
moments Ci = (n*)/(n)* can be expressed as an expan¬ 
sion in the inverse number of average participants l/(fVp), 


a = 


Np 

(Np) 


Nc , 


(1 


o 



(14) 


where the average ratio of hard to soft component 


{Nc){nc) 

{Np) (rip) 


(15) 


depends on the elementary nucleon-nucleon dynamics 
and the nuclear geometry. If we do not consider periph¬ 
eral interactions alone, {Np),{Nc) >> 1, we can then 
consider only the zeroth order in the expansion (14). In 


2 








this case, only the average ratio of hard to soft compo¬ 
nent remains, the other dynamics of elementary soft and 
hard processes hidden in {n^) and (n*) with i > 2 is 
washed away by the nuclear geometry. 

When the hard contribution can be neglected, namely 
a; —*■ 0, the normalized multiplicity moments are just the 
normalized participant moments, 

(TV*) 

O = C., = ^ . (16) 

This is the case discussed in Ref. at SPS energies. 


III. NUCLEAR GEOMETRY AND ENERGY 
DEPENDENCE OF HARD CONTRIBUTION 

let’s first determine the soft and hard components (rip) 
and (ric) in elementary nucleon-nucleon collisions. To 
this end, we compare the average multiplicity with the 
experimental data for central Au—Au collisions. Since we 
did not introduce rapidity dependence in our discussion, 
we consider only the central rapidity region where the 
data show a plateau structure for different centrality bins. 
By comparing the average participant number (fVp), the 
average multiplicity per participant pair 

fh!|j = 2{V(l+x|, (17) 

and the multiplicity for PP 

{npp) = 2{np) + {ric} (18) 

with the experimental data |^,|^ in the central rapidity 
region \t]\ < 1 at RHIC and the parametrization of the 
PP data Q 

(n)pp = 2.5-0.25 In s-b 0.02310^5 , (19) 

we can determine at different energies the average ra¬ 
tio X and the minimum participant number Nmin which 
is used to select centrality in calculating geometry mo¬ 
ments. Using a Wood-Saxon distribution, 

PAir) = -, [ (firpAir) = A , (20) 

1 -b e o 



(Np) 

(n) 

0.b{NA 

^min 

X 

56 

330 

2.47 

287 

0.32 

130 

343 

3.24 

312 

0.58 

200 

344 

3.78 

314 

0.74 


TABLE I. The ratio x of hard to soft component and the 
geometry parameter Nmin determined from the eomparison 
with the data of eentral An — Au collisions at RHIC. 


with the parameters a = 0.53/to, Ra = I.IA^^^fm for 
^^'^Au and taking = 37 mb at ^/s = 56 A GeV (ctat = 
41 mb for -/s = 130, 200 A GeV) Q, the two parameters 
are shown in Tab.(^. We see that at RHIG energies, 
a: < 1, the soft component is still more important than 
the hard component. 

The influence of nuclear geometry is twofold: The av¬ 
erage numbers (Np) and (Nc) and the fluctuations of Np 
and Nc around their average values. For central colli¬ 
sions the average numbers {Np) and (Nc) are huge, but 
the fluctuations are small. This can be seen clearly in 
Tab.(|) where (Np) > 330 and 287 < Np < Np{b = 0). 
For minimum-bias events the average numbers are rela¬ 
tively small, but the fluctuations are the maximum. 

The multiplicity (n) is only related to the average val¬ 
ues (Np) and (A/). When the hard contribution van¬ 
ishes, the average multiplicity is proportional to {Np). 
The hard contribution reflected in the ratio x leads to 
an extra {Nc) dependence. The centrality dependence of 
the average multiplicity per participant pair ( 0 ) can be 
calculated by changing the minimum participant number 
Nmin from 2 to Np{b = 0). In Fig.(^ it is compared with 
the data in the central rapidity region [ryj < 1 for the cen¬ 
tral Au — Au collisions at y/s = 130 A GeV [^. The 
extra geometry dependence induced by the hard compo¬ 
nent is weak. 
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FIG. 1. The eentrality dependence of the average multiplic¬ 
ity normalized to per participant pair and its comparison with 
the RHIC data. 

Since {N;) = {{Np/{Np)Y){Np)^ and (TV*) = 

{{Nc/{Nc))^){Nc)'^, the multiplicity moments (n®) for 
i > 2 are associated with both the average numbers 
{Np) and {Nc) and the fluctuations in Np and Nc. From 
Eq.(|l^ the normalized moments Ci depend on the fluc¬ 
tuations and the average ratio x of hard to soft compo¬ 
nent. Fig.(^ shows the centrality and energy dependence 
of a;. At any energy the centrality dependence is very 
weak. Therefore, the behavior of the normalized mo¬ 
ments Ci is mainly controlled by the fluctuations in Np 
and Nc. Let’s first consider the limit of no fluctuations. 


(dN/di7)/((Np>/2) 



Vs=130GeV, |77l<1 
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Np = (Np), Nc = (Nc). In this limit, 


p{Np) = 6np{n^) , (21) 

we have 

Q = Qp = 1 . (22) 

In this case there is no difference between the two- 
component and one-component model. Although fluctu¬ 
ations around the average numbers always exist, and it is 
difficult to choose events with the same impact parameter 
b, namely with the same Np and Nc, in experiments, for 
very central collisions with large {Np) and {Nc), Np and 
Nc fluctuate in a narrow region, the case is then similar 
to the above limit. 



<Np) 

FIG. 2. The energy and centrality dependence of the aver¬ 
age ratio X of hard to soft component. 

The fluctuations grow up when the minimum partic¬ 
ipant number Nmin decreases from its maximum value 
Np{0). Fig.(^) shows the centrality dependence of the 
fluctuations {NpNi)/{Np)''{Nc)^■ As the orders i and j 
are not too small, the fluctuations are very strong for 
minimum-bias events. 
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FIG. 3. The centrality dependence of the geometry fluctua¬ 
tions. 


In order to see the contribution from the hard pro¬ 
cesses, we define the ratio of the normalized moments 
with and without consideration of the hard component. 


n = 



(23) 


The centrality and energy dependence of is shown in 
Fig.(^). While there is no remarkable difference between 
Cip and Ci in central collisions, the big fluctuations in Np 
and Nc in minimum-bias events enhance the hard contri¬ 
bution, and this enhancement become more and more im¬ 
portant when colliding energy increases. At -y/s = 200 A 
GeV, the hard contribution to C 5 is almost 50%. 
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FIG. 4. The ratio of the two-component to one-component 
normalized moment as a function of the centrality. 


IV. COMPARISON WITH EFFECTIVE MODEL 
WITHOUT EXPLICIT HARD COMPONENT 

The effect of the hard scattering processes on the av¬ 
erage multiplicity can be effectively described in the one- 
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component model by modifying the participant number 



, Np^N;^ , a>l . (24) 

By comparing the average multiplicity (n) = {Np){np^^) 
with the RHIC data listed in Tab.(|), we can determine 
the power a and the average contribution of each effective 
soft source Corresponding to the colliding energy 

yi = 56,130,200 A GeV, we have a = 1.04,1.07,1.09, 
respectively. 

In the effective one-component model, the normalized 
moments are just the effective participant moments, 


= 


{Np 

{N^Y 


(25) 
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FIG. 5. The ratio of two-component to ef¬ 
fective one-component normalized moment as a function of 
the centrality. 


when the peripheral interactions are not considered 
alone. While the contribution of the hard processes to the 
average multiplicity through the average binary collision 
number {NY} can be equivalently expressed by increasing 
the average participant number from {Np) to {N^), the 
fluctuations in Nc can not be effectively included in the 
fluctuations in (A^“). This can be seen clearly in Fig.(||) 
which shows the ratio 


Ri = 



(26) 


as a function of the centrality for Au — Au collisions. 
From the comparison with Fig.(||), Ri < r^, the fluctua¬ 
tions in Nc are partly included in the fluctuations in the 
effective participant number A^“. However, the differ¬ 
ence between the two-component model and the effective 
one-component model is still remarkable in minimum- 
bias events, especially for the higher order moments and 
at high energies. 


V. CONCLUSIONS 

The huge average participant number {Np) and bi¬ 
nary collision number {Nc) in relativistic heavy ion col¬ 
lisions make it difficult to extract dynamic information 
on hard processes from the geometry background. Dif¬ 
ferent from the multiplicity moments (n*) which depend 
on both the average numbers {Np) and {Nc) and the fluc¬ 
tuations in Np and Nc strongly, the normalized moments 
Ci = {'N)/{nY have only weak {Np) and {Nc) depen¬ 
dence, and are mainly associated with the fluctuations in 
Np and Nc- Therefore, the geometry background for Ci 
is not so complicated as that for (n*). 

We have investigated the normalized moments Ct in 
the frame of a general two-component model. When the 
hard component can be neglected at SPS energies, Ct are 
completely determined by the geometry fluctuations, the 
dynamics is totally washed away. When the hard pro¬ 
cesses become important at RHIC energies, the average 
ratio of hard to soft component depends on the centrality 
weakly, and Ci are dominated by the fluctuations. For 
central collisions where the fluctuations are weak, Ci ap¬ 
proach to 1, the dynamic information can not be seen 
in Ci- However, the big fluctuations in minimum-bias 
events make us to see clearly the difference between the 
models with and without hard component. 

While the average effect of the hard processes can be ef¬ 
fectively described in the one-component model by mod¬ 
ifying the participant number, we have found that the 
fluctuations in the binary collision number can not be 
fully included in the fluctuations in the effective partici¬ 
pant number. 
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